Optimal routing on complex networks 
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We present a novel heuristic algorithm for routing optimization on complex networks. Previously 
proposed routing optimization algorithms aim at avoiding or reducing link overload. Our algorithm 
balances traffic on a network by minimizing the maximum node betweenness with as little path 
lengthening as possible, thus being useful in cases when networks are jamming due to queuing 
overload. By using the resulting routing table, a network can sustain significantly higher traffic 
without jamming than in the case of traditional shortest path routing. 
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The ever-increasing amount of information and goods 
transported along complex networks raises the question 
about their ultimate transport capacity. Generally, the 
transport capacity of a network is limited by two factors: 
link capacity (bandwidth in the case of the Internet or 
wireless networks) and node (router) latency. Tradition- 
ally, both fixed and ad-hoc network routing is based on 
the idea of maintaining a table of the shortest paths (or 
the best available approximation to the shortest paths) 
between any two nodes of the network and forwarding 
the information packages along these paths Q, ■ The 
length of a path is computed as the sum of the weights 
assigned to the links that form the path. In the case 
of the Internet, link weights are typically assigned man- 
ually by operators according to simple rules based on 
experience Q . Recently, a series of heuristic algorithms 
have also been proposed for network traffic optimization 
|H H H 11 0] These rules and algorithms are aimed 
at avoiding or reducing link overload by a judicious link 
weight assignment. The cost of each link is assessed as 
a monotonically increasing function of the ratio between 
traffic and capacity and then the weights are adjusted to 
minimize the sum of the costs of all links. 

This approach, however, does not take into account 
delays due to node latency (information package pro- 
cessing time) and has the disadvantage that too many 
of the shortest paths pass through a few nodes, called 
hubs. As a result, in situations of high network traf- 
fic, these hubs will experience congestion (long message 
queues) and eventually jamming, causing the network to 
break apart in a multitude of disconnected subnetworks. 
In light of this behavior, the optimality of the shortest 
path routingas currently implemented has been recently 
questioned 0, @, H El Hi El El Q E3- It has 
been shown, for example, that dynamic routing proto- 
cols which allow for a certain degree of stochasticity or 
take into account the congestion status of the nearest 
neighbors significantly improve the transport capacity of 
a network El EH El El D El • A more systematic 
approach is to find better static (strictly table-driven) 
routing protocols that avoid the hubs whenever possible 



and convenient (i.e. when avoiding a hub does not lead 
to congestion on another node). Recent studies H 
have shown that this is possible and propose new rout- 
ing algorithms which lead to improved transport capacity 
(quantified by the packet insertion rate at which jamming 
occurs). An open question is how much larger the actual 
transport capacity can be than the results presented in 
Refs. HQ- -*- n tms Letter, we show that significant im- 
provement in the transport capacity of a network can be 
achieved by systematically adjusting the traffic routing 
to minimize the maximum betweenness on the network. 
Our algorithm leads to higher transport capacity than 
those presented in Refs. H 0|- The transport capac- 
ity also exceeds the analytical estimate for its maximum 
value given in [jj. Furthermore, we argue that our al- 
gorithm achieves near-optimal routing for uncorrelated 
scale-free networks. 

To facilitate comparison, we use the same network 
model as in Ref. Q- All results presented in this Letter 
are for undirected, uncorrelated scale-free networks with 
an exponent of the power-law degree distribution 7 = 2.5, 
generated using the configuration model. The number of 
nodes N varies between 25 and 1600. For simplicity, we 
assume that all nodes have the same processing capac- 
ity of 1 package per time step and that new information 
packages are inserted at every node at the same average 
rate of r packages per time step. The destinations of the 
packages are chosen at random from among the other 
N — 1 nodes on the network. However, the algorithm can 
be generalized for nodes with different processing capac- 
ities and for arbitrary traffic demands. Given a routing 
table, the betweenness Bi of node i is defined as the 
sum of all fractional paths that pass through that node. 
The fraction of times a message passes through node i on 
its way from a source node s to a target node t is com- 
puted as follows: the source node s is assigned a weight 
1 and then the weight of every node along each path is 
split evenly among its predecessors in the routing table 
on the way from f to s and added to the weights of the 
predecessors. The average number of packages passing 
through a given node i is then < w >;= rBi/(N — 1). 
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FIG. 1: (Color online) Ensemble average of the network max- 
imum betweenness as a function of the number of nodes for 
four routing protocols. 

Jamming occurs at the critical average insertion rate r c 
at which the average number of packages processed by 
the busiest node reaches unity. Consequently, r c is given 
by 



where B max is the highest betweenness of a node on the 
network. Thus, to achieve optimal routing, the highest 
betweenness B max should be minimized. An important 
point is that, even though the minimization procedure 
pertains to a single scalar quantity, such an optimization 
algorithm will implicitly reshape the betweenness land- 
scape across the whole network, lowering traffic through 
the initially busy nodes at the expense of increased traffic 
through the initially idle nodes until the traffic spreads 
out and an as narrow as possible betweenness distribu- 
tion is achieved. 

The problem of finding the exact optimal routing is 
mathematically tied to the problem of finding the mini- 
mal sparsity vertex separator 0, which has been shown 
[l7| to be an iVP-hard problem. This means that the 
number of flops necessary for the computation of an ex- 
act solution increases with the number of nodes N faster 
than any polynomial. We propose a heuristic algorithm 
which finds near-optimal solutions for the routing prob- 
lem in time 0(N 4 ) at worst (0(N 3 ) for one iteration 
and requiring O(N) iterations), fn its simplest form, the 
algorithm proceeds as follows: 

1. Assign weight I to every link and compute the short- 
est paths between all pairs of nodes. 

2. Compute the betweenness of every node. 

3. Find the node which has the highest betweenness 
Bmax and add I to the weight of every link that connects 
it to other nodes. 

4. Recompute the shortest paths. Go back to step 2. 
Note that the algorithm picks the "least fit" element of 
a set and changes its parameters. Therefore, it is a from 
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FIG. 2: (Color online) Ensemble average of the network av- 
erage betweenness as a function of the number of nodes for 
four routing protocols. 

extremal optimization [l^. However, this algorithm as- 
signs parameters in a deterministic way, unlike many of 
the other existing extremal optimization algorithms. 

Before presenting the results, we note that the average 
betweenness B avg on a given network using a given rout- 
ing table provides an absolute lower bound for the max- 
imum betweenness. This lower bound would be achiev- 
able only if one could optimize routing to the point where 
all nodes experience the same traffic. The difference be- 
tween the maximum and average betweenness is also a 
measure of the width of the betweenness distribution. 
Moreover, the average betweenness in the case of short- 
est path routing with all weights set to I (which in the 
remainder of this Letter will be called shortest path rout- 
ing) constitutes the lower bound for the average between- 
ness computed using any arbitrary routing table, since 
any changes from the shortest path routing (including 
those resulting from an optimization algorithm) will re- 
sult in longer paths, thus adding to the sum of all be- 
tweennesses on the network. It is thus apparent that a 
good optimization algorithm is required to have at least 
two properties: (I) minimize the difference between the 
maximum and average betweenness, and (2) do this while 
keeping the difference between the average betweenness 
computed using the optimized routing table and the one 
computed using the shortest path routing table as low as 
possible. 

In the following, networks of a given size N are char- 
acterized by the ensemble averages of the maximum be- 
tweenness < B max > and average betweenness < B avg > 
computed over a set of 100 realizations of the network. 
Computer simulation results presented in Ref. |7f suggest 
a power-law functional dependence of < B max > on the 
network size N. Our results confirm this power-law de- 
pendence for < B max > and show a similar dependence 
in the case of < B avg >. 

Results for the ensemble average of the maximum be- 
tweenness < B max > as a function of the network size N 
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FIG. 3: Maximum betweenness as a function of the number 
of iterations for a network with 196 nodes. 
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FIG. 4: (Color online) Distribution of node betweennesses be- 
fore (black shaded bins in (a) and black crosses in (b)) and 
after (red hollow bins in (a) and red circles in (b)) optimiza- 
tion for a network with 400 nodes. 



for four routing protocols are presented in Fig. 1 . Similar 
results for the ensemble average of the network average 
betweenness < B avg > are shown in Fig. 2. The results 
for the shortest path (SP) routing are represented by cir- 
cles, while those obtained for the routing given by our 
optimization algorithm (OR) are represented by squares. 
For comparison, we also show results obtained using the 
efficient routing (ER) and hub avoidance (HA) protocols 
described in Refs. [(j and Q, represented by diamonds 
and triangles respectively. The exponents resulting from 
fitting the data points in each set are given in Table 1, 
where the quoted errors are 2er estimates. 

It is apparent from Figs. 1 and 2 that the routing based 
on our optimization algorithm leads to a far smaller dif- 
ference between the maximum and the average between- 



ness than in the case of the other three routing proto- 
cols. Moreover, the maximum and the average between- 
ness scale with N with the same exponent, which is a 
strong argument in favor of the optimality of the rout- 
ing. Finally, the difference < B avg > r - < B avg >sp, 
while larger than < B avg >er — < B avg >sp or 

< B avg >ha - < B avg > S p (which is explained by 
the need to have slightly longer paths around the hubs), 
is kept quite low and < B aV g >or scales with the net- 
work size ./V with an exponent only slightly higher than 

< B avg >sp. 

As expected with a heuristic algorithm, the evolution 
of the maximum betweenness as a function of the number 
of iterations is not monotonic, but exhibits a decreasing 
trend and eventually the maximum betweenness "con- 
verges" in the sense that it becomes confined to a narrow 
band. This is exemplified in Fig. 3, which is a plot of 
Bmax versus the number of iterations for a given network 
with 196 nodes. The algorithm can be implemented to 
pick the configuration corresponding to the global mini- 
mum of B max . 

Fig. 4 provides insight into how the algorithm works by 
comparing the initial and final betweenness distributions 
in the case of a network with 400 nodes. Fig. 4a shows 
histograms of the betweenness distribution before and af- 
ter optimization, while Fig. 4b shows the betweennesses 
plotted against vertex index. Initially, the majority of the 
nodes have very low betweenness, but a small number of 
them are spread over a very wide betweenness range. Af- 
ter optimization, all node betweennesses are confined to 
a narrow band, whose upper edge is quite well defined. 
Most of them are uniformly distributed within the band, 
but there is a very sharp peak at the upper edge. There is 
a significant decrease in the number of very low between- 
ness nodes. A plot of the final (OR) betweenness versus 
the initial (SP) betweenness in the case of a network with 
400 nodes is shown in Fig. 5. It is apparent from this plot 
that the algorithm performs remarkably well, by lower- 
ing the traffic through all nodes whose initial (shortest 
path) betweenness lies above a certain critical value until 
they all reach essentially the same critical betweenness. 
On the other hand, virtually all nodes whose initial be- 
tweenness lies below the critical value experience higher 
traffic, many of them (especially those with higher ini- 
tial betweenness) reaching the critical value. It is still 
an open question whether an improved algorithm can 
achieve a lower critical betweenness by further raising 
the traffic through some initially low betweenness nodes. 
On the other hand, it is clear that not all low between- 
ness nodes can have their betweenness increased without 
unduly lengthening paths or increasing traffic through 
other nodes which are prone to congestion. The simplest 
examples are those of a small subnetwork whose only 
connection to the rest of the network is through a single 
link to a high degree (or otherwise high SP betweenness) 
node, or a triangle connected to the rest of the network 
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1.088 ±0.019 
1.634 ± 0.010 


1.186 ±0.024 
1.185 ±0.009 


1.165 ±0.009 
1.315 ±0.017 


1.080 ± 0.009 
1.542 ±0.010 



TABLE I: Exponents of the < B avg > and < B max > power-law scaling with network size N. 



work to bear significantly higher traffic than in the case 
of shortest path routing. Network capacity is improved 
by a factor which increases with network size according 
to a power law. 

The authors acknowledge support from the NSF 
through grant No. DMR-0427538 and also from SI Inter- 
national through A. Williams of the Air Force Research 
Laboratory Information Directorate under contract No. 
FA8750-04-C-0258. 



FIG. 5: Correlation plot of the final (OR) versus initial (SP) 
betweenness for a network with 400 nodes. 



only by containing such a node. In these cases, there is no 
way of diverting any of the high SP betweenness node's 
traffic between other nodes through the structures men- 
tioned above. The latter will have low betweenness even 
in the case of rigorously optimal routing. 

The difference between the SP and OR distribution of 
the travel times between the various nodes in situations 
when the SP routing does not lead to jamming is still an 
open question and will be the subject of a future study. 
However, we argue that, at least in situations of con- 
gested traffic, most travel times are shorter in the case 
of OR routing. It is known from queuing theory that 
the average queue length < q > is given (assuming unity 
processing power) by 0, El 



< q >- 



<w> 
1- < w > : 



(2) 



with < w > defined above Eq. 1. Due to this strongly 
nonlinear relationship, which diverges as < w > ap- 
proaches unity, it seems reasonable to assume that by 
avoiding the passage through hubs with very high be- 
tweenness most travel times become shorter, in spite of 
the fact that the routes pass through more nodes. This 
conclusion is also supported by the results for the average 
path length and average travel time presented in Ref. [6j . 

In summary, we have presented a simple heuristic al- 
gorithm for routing optimization on complex networks 
and demonstrated its usefulness for scale-free networks. 
This algorithm is useful in situations when network jam- 
ming is primarily due to queuing overload. Our results 
show that the application of this algorithm allows a net- 
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